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Abstract. We study gradient estimates of g-harmonic functions u of the frac- 
tional Schrodinger operator A Q / 2 + q, a £ (0, 1] in bounded domains D C K, d . 
For nonnegative u we show that if q is Holder continuous of order r\ > 1 — a then 
Vu(x) exists for any x £ D and |Vu(x)| < cu(x)/(dist(x,dD) A 1). The exponent 
1 — a is critical i.e. when q is only 1 — a Holder continuous Vzt(s) may not exist. 
The above gradient estimates are well known for a £ (1,2] under the assumption 
that q belongs to the Kato class 2 a ~ 1 ■ The case a £ (0, 1] is different. To obtain 
results for a £ (0, 1] we use probabilistic methods. As a corollary, we obtain for 
a £ (0, 1) that a weak solution of A Q / 2 u + qu = is in fact a strong solution. 

CO' 



1. Introduction 

Let a G (0,2), del and q belong to the Kato class 3 a . We say that a Borel 
function u on R d is q-harmonic in an open set D C R d iff 



9h 



m(x) = 



exp ^ g(X s ) rfs^j m(X Tw ) 



x e W, (l) 



for every open bounded set W, with W G D. Here X 4 is the symmetric a-stable 
process in R d , tw the first exit time of X t from W, and we understand that the 
expectation in (JT} is absolutely convergent. 

It is possible to express the above probabilistic definition in analytic terms. Namely, 
it is known [SI Theorem 5.5] that if u is g-harmonic in open set D C R^ then u is a 
weak solution of 

A a/2 u + qu = 0, on D. (2) 

Here A a//2 := — (— A) Q//2 is the fractional Laplacian. On the other hand if D C R d 
is an open bounded set and (D, g) is gaugeable then a weak solution of d2J is a 
g-harmonic function on D after a modification on a set of Lebesgue measure zero 
(for more details see Preliminaries). 

It is known [8] that if u is g-harmonic in D then it is continuous in D. The 
purpose of this paper is to derive further regularity results of g-harmonic functions. 
The main result is the following. 

Theorem 1.1. Let a G (0, 1], d G IN and D C R d be an open bounded set. Assume 
that q : D — > R is Holder continuous with Holder exponent rj > 1 — a. Let u be 
q-harmonic in D. If u is nonnegative in R d then Vu(x) exists for any x G D and 
we have 

u( x) 

\Vu(x)\<c - y x E D, (3) 
5 D (x) A 1 

where Su(x) = dist(x, 3D) and c = c(a, d, rj, q). 
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If u is not nonnegative in R d but ||w||oo < 00 ^en V«(x) exists for any x G D 
and we have 



Vu{x)\<c 1 °° x e D, (4) 

dn X A 1 




where c = c(a, d, 77, g) . 

The existence of V«(i) and similar gradient estimates are well known in the 
classical case for a = 2, see e.g. [16] and for a G (1,2), see [10j. These results for 
a G (1,2] were shown under the assumption that q G 2 a ~ 1 ■ The biggest difference 
between the cases a G (0,1] and a G (1,2] is the fact that for a G (0,1] the 
function y — > \V x Gd(x, y)\ is not integrable while for a G (1,2] is integrable. Here 
Gd(x,v) is the Green function for A"/ 2 with Dirichlet condition on D c . The fact 
that y — > \V x Gd{x, y)\ is integrable was widely used in [10J for a G (1,2), see e.g. 
[TO] Lemma 5.2]. For a G (0,1] more complicated method must be used. Key 
ingredients of the method for a G (0, 1] may be briefly described as the combination 
of some estimates of the Green function and some self-improving estimates used 
in the proof of Theorem 11.11 The proof of the estimates of the Green function 
is mainly probabilistic. It is based on the representation of symmetric a-stable 
processes as subordinated Brownian motions and the reflection principle for the 
Brownian motion. This probabilistic idea is similar to the one used in the paper 
by B. Bottcher, R. Schilling, J. Wang, where they study couplings of subordinated 
Brownian motions, see Section 2 in [TT]. More remarks about these probabilistic 
methods are at the end of Section 3. 

From analytic point of view Theorem 11.11 gives some regularity results for weak 
solutions of (f2]). It is worth to notice that regularity results of weak solutions 
of equations involving the fractional Laplacian have attracted a lot of attention 



One may ask whether it is possible to weaken the assumption in Theorem 1 1 . 1 1 1 hat 
q is Holder continuous with Holder exponent rj > 1 — a. It occurs that the exponent 
77 = 1 — a is critical in the following sense. 

Proposition 1.2. For any a G (0, 1], d G IN and any open bounded set D C R rf there 
exists q : D — > [0, 00) which is 1 — a Holder continuous, a function u : R d — > [0, 00) 
which is q-harmonic in D and a point z G D such that Vu(z) does not exist. 

The proof of this proposition is based on the estimates of the Green function of 
the killed Brownian motion subordinated by the a/2-stable subordinator. These 
estimates were obtained by R. Song in [27] . 




When a g-harmonic function u vanishes continuously near some part of the bound- 
ary of D and D C R d is a bounded Lipschitz domain then the estimates obtained 
in Theorem 11.11 are sharp near that part of the boundary. 

Theorem 1.3. Let a G (0, 1], d G W ; D C R d be a bounded Lipschitz domain and 
q : D — >■ R be Holder continuous with Holder exponent i] > 1— a. Let V C R d be open 
and let K be a compact subset ofV. Then there exist constants c = c(D, V, K, a, q, rf) 
and e = e(D, V, K, a, q, r/) such that for every function u : H d — >■ [0, 00) which is 
bounded on V , q-harmonic in D C\V and vanishes in D c n V we have 



recently, see e.g. [17] . [26] . 




-1 



u(x) 



< |Vtt(x)| < c 



u(x) 



x G K n D, 8 D (x) < e. 



c 



S D (x) 
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Similar result was obtained for a = 2 in [4 J and for a G (1, 2) in [10], see Theorem 
5.1. 

As an application of our main result we obtain gradient estimates of eigenfunc- 
tions of the eigenvalue problem of the fractional Schrodinger operator with Dirichlet 
boundary conditions. These estimates are formulated and proved in Section 6. 

As another application of our main result we show for a G (0, 1) that under some 
assumptions on q a weak solution of A a ^ 2 u + qu = is in fact a strong solution. Note 
that in the following corollary we do not have to assume that (D, q) is gaugeable. 

Corollary 1.4. Let a G (0, 1), d G IN and D C R d be an open bounded set. Assume 
that q : D — )■ R is Holder continuous with Holder exponent 77 > 1 — a and either 
u is nonnegative on R d or \\u\\oo < 00. If u is a weak solution of |U) then (after a 
modification on a set of Lebesgue measure zero) u is continuous on D and it is a 
strong solution of (TJ)]. 

The paper is organized as follows. Section 2 is preliminary; we collect here basic 
facts concerning the fractional Laplacian, the fractional Schrodinger operator and 
g-harmonic functions. In Section 3 using probabilistic methods we obtain estimates 
of the Green function, which will be essential in the rest of the paper. In Section 4 
the main result of the paper is proved. Section 5 contains proofs of Proposition 11.21 
and Theorem 11.31 Section 6 concerns applications of the main result. 

2. Preliminaries 

Most of the terminology and facts presented here are taken from [S] and [9]. The 
notation c(a, b, . . .) means that c is a constant depending only on a, b, . . .. Constants 
are always positive and finite. We adopt the convention that constants may change 
their value from one use to another. As usual we write x Ay = min(x, y), x V y — 
m&x(x,y) for x,y G R, | )| for any function u : R d — > R, 

B(x, r) = {y G R d : \x — y\ < r} for x G R d , r > 0. By e,, i — 1, . . . , d we denote 
the standard basis in R d . 

We denote by (X t ,P x ) the standard rotation invariant ("symmetric") a-stable 
process in R d , a G (0, 2] with the characteristic function E° exp(i£X t ) = exp(— £|£| a ), 
£ G R d , t > 0. E x denotes the expectation with respect to the distribition P x of 
the process starting from x G R d . We have P x (X t G A) — j A p(t,x,y) dy, where 
p(t, x, y) = pt{y — x) is the transition density of X t . 

For a < d the process X t is transient and the potential kernel of X t is given by 

K a (y-x)= ! p(t,x,y)dt= x,yeR d , (5) 

Jo \y ~ x \ 

where A(d,j) = T{{d - ^)/2)/{2^n d / 2 \r{-f /2)|) [5J. When a > d the process is re- 
current and it is appropriate to consider the so-called compensated kernels. Namely 
for a > d we put 

/■oo 

K a {y-x)= \ (p(t,x,y)-p(t,0,x ))dt, 
Jo 

where Xq — for a > d — 1, xq — 1 for a — d — 1 and xq = (0, 1) for a — d — 2. 
For a. — d — 1 we have 

K a (x) = - log 



7T \\X 
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For any open set D C R d we put td = inf{t > : X t ^ D} the first exit time 
of X t from D and we denote by pn{t,x,y) the transition density of the process X t 
killed on exiting D. The transition density is given by the formula 

p D (t,x,y) = p(t,x,y) - E x (p(t- T D ,X(T D ),y),r D <t), x,y G D, t > 0. 

We put poit, x, y) = if x G D c or y G -D c . It is known that for each fixed t > the 
function "j ') is bounded and continuous on D x D. When d > a and D G R d 
is an open set or d = 1 < a and D C R d is an open bounded set we put 



poo 

G D (x,y) = p D (t,x,y)dt, x,y G D, 
Jo 



Gd(x, y) = if x G D c or y G -D c . We call GdI^, 2/) Green function for D. It is 
known that Gd(x, •) is continuous onD\ {x}. For any open bounded set D C H d 
we define the Green operator Gd for D by 



G D f(x) = j G D (x,y)f(y)dy. 
We assume here that / is a bounded Borel function / : D — > R. We have 

G D f(x)=E x r f(X s )ds. 



JO 

Now we briefly present basic definitions and facts concerning the fractional Lapla- 
cian and the fractional Schrodinger operator. We follow the approach from [8]. We 
denote by L 1 the space of all Borel functions / on R d satisfying 

_\m\ 

For / G L 1 and x G R d we define 



/ 



. , dx < DC. 
1 + \x\) d+a 



A a ^ 2 f(x) = A(d, -a) lim [ f ^ dy, 

K J 1 ' ^J\y-x\>e \y~x\ d +« 

whenever the limit exists. 

We say that a Borel function q : H d — > R belongs to the Kato class 3 a iff q satisfies 



lim sup / \q(y)K a (y- x)\ dy = 0. 

r -l" xeR d J\y-x\<r 

For any a G (0, 2), q G 3 a we call A Q//2 + g t/ie fractional Schrodinger operator. 

Let a G (0,2), q E d a and D C R d be an open set. For u G L 1 such that 
«g G Ll oc (D) we define the distribution (A Q / 2 + g)w in D by the formula 

((A a/2 + q)u, tp) = (it, A^V + y € Cf(D), 

(cf. Definition 3.14 in [8]). We will say that u is a weak solution of 

(A a / 2 + q)u = (6) 

on D iff it G uq G L\ 0C (D) and (jSJ) holds in the sense of distributions in D. We 
will say that u is a strong solution of © on D iff w G £ 1 , uq G L\ oc (D) and (J6]) holds 
for any x & D. 

For a G (0,2), g G £P the multiplicative functional e 9 (i) is defined by e 9 (£) = 
exp f J * g(X s ) cfej , t > 0. For any open bounded set D C R d the function 

u D {x) = E x (e q (r D )) 
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is called the gauge function for (D, g); when it is bounded in D we say that (D, g) 
is gaugeable. There are several other equivalent conditions for gaugeability, in par- 
ticular there is a condition in terms of the first Dirichlet eigenvalue of A Q//2 + q on 
D, see below. 

Let u be a Borel function on R d and let q E 3 a - We say that u is q-harmonic in 
an open set D C R d iff 

^I^'Mt^J], xgW, (7) 

for every bounded open set W with W C D. u is called regular q-harmonic in Z) iff 

= [e ? (rD)w(X TD ); r D < oo] , x G D. (8) 

We understand that the expectation in (0) and (JSJ) is absolutely convergent. 

By the strong Markov property any regular g-harmonic function in D is a in- 
harmonic function in D. By [HI Theorem 4.1] any g-harmonic function in D is 
continuous in D. By [91 (4.7)] any g-harmonic function in D belongs to L 1 (when 
D 7^ 0). It follows that if u is a g-harmonic function in D then uq G L\ oc (D). 

Let a 6 (0,2), g 6 £J a . If w is a g-harmonic function in an open set Dctf then 
it is a weak solution of (A Q//2 + q)u = on D. Conversely assume that D C R d is 
an open bounded set and (D, g) is gaugeable. If a function u is a weak solution of 
(A Q / 2 + g)w = 0onf then after a modification on a set of Lebesgue measure zero, 
u is g-harmonic in D (see [HI Theorem 5.5]). 

It is known that if u is g-harmonic in open set D C R d then, unless w = on 
D and w = a.e. on D c , (W, q) is gaugeable for any open bounded set W such 
that W C D (see [HI Lemma 4.3]). We will often use the following representation 
of g-harmonic functions. If u is g-harmonic in an open set D C JR d then for every 
open bounded W with the exterior cone property such that W C D we have 

u(x) = E x u(X TW ) + G w (qu)(x), xeD. (9) 

This follows from [91 Proposition 6.1] and continuity of g-harmonic functions. 

By saying that q : D — > R, is Holder continuous with Holder exponent rj > 
we understand that there exists a constant c such that for all x, y G D we have 
|g(x) — g(y)| < c|x — y\ v . 

We finish this section with some basic information about the spectral problem for 
A a / 2 + g. Assume that D C R d is an open bounded set, a G (0,2), g G d a - Let us 
consider the eigenvalue problem for the fractional Schrodinger operator on D with 
zero exterior condition 

A Q/ V + g<^ = -\<p on D, (10) 

ip = on D c . (11) 

It is well known that for the problem ( TTQTfTTj) there exists a sequence of eigenvalues 
{A n }~ =1 satisfying 

Ai < A2 < A3 < . . . , lim A n = 00, 

and a sequence of corresponding eigenfunctions {</?n}5£Li, which can be chosen so 
that they form an orthonormal basis in L 2 (D). All tp n are bounded and continuous 
on D and tp\ is strictly positive on D. It is also well known that gaugeability of 
(D,q) is equivalent to Ai > see [12] Theorem 3.11], cf. [151 Theorem 4.19]. We 
understand that f llip holds for all x G D c and (1101) holds for almost all x G D. The 
eigenvalue problem ( TTUlfTT]) was studied in e.g. [T2], [22] and very recently in [2D] . 
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For more systematic presentation of the potential theory of fractional Schrodinger 
operators we refer the reader to [5] or to [7j. 

3. Estimates of the Green function 
In this section we fix i G {1, . . . , d} and use the following notation 

H = {(y l} ...,y d )eR d : Vi >0}, 

H = {(y 1} ...,y d )eR d : y l = 0}. 
Let R : R d — > R d be the reflection with respect to H . For any x G R d we put 

x = R(x). 

We have x = x — 2xjej, where (ei, . . . , e^) is the standard basis in R d and x = 
(xi, . . . , Xd)- We say that a set D C R d is symmetric with respect to Hq iff R{D) = D. 
For any set D C R d , which is symmetric with respect to Hq we put 

D + = {(y u ...,y d )eD: Vl >0}, D_ = {(y u . . . , y d ) G D : y< < 0}. (12) 

Let -Bf be the d- dimensional Brownian motion starting from x G R d (with the 
transition density (4iit)~ d / 2 e~\ x ~ y \ /( 4 *)) and ^ be the a/2-stable subordinator start- 
ing from zero, a G (0, 2), independent of B t (E~ s ' nt = e~ tsa/2 ). It is well known that 
the (i-dimensional symmetric a-stable process X t , a G (0,2), starting from x G R d 
has the following representation 

X t = B Vt . 

Let 

T = inf{s > : B s G H }. 
Assume that the Brownian motion B t starts from x G H. We define 

R(B t ) for t < T 



B t 



B t for t > T. 



That is B t is the mirror reflection of B t with respect to Hq before T and coincides 
with B t afterwards. It is well known that B t is the Brownian motion starting from 
x. Now set 

x t = B Vt . 

X t is the symmetric a-stable process starting from x. The above construction is 
taken from Section 2 in [11]. When discussing probabilities of B t , B t , i] t , X t , X t we 
will use P§, P§, P v , P x , P x respectively. 

Now we need to consider another process, which is a subordinated killed Brownian 
motion. We define it as follows 

X t = (B H )ri t , 

where Bf is the Brownian motion B t (starting from x G H) killed on exiting H 
and i] t is the a/2-stable subordinator starting from zero, independent of B t . When 
discussing probabilities of Bf , X t we will use P^h, P x respectively. The general 
theory of subordinated killed Brownian motions was studied in [28J . 

For any open set D C R d , which is symmetric with respect to H we put 

f D+ = inf{t > : X t i D + }, 

where D + is given by f)12p . 
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By Pd + (t, x, y) we denote the transition density of the process X t killed on exiting 
D + . The idea of considering pD + (t, x, y) comes from [U Section 4]. 

Recall that pD(t,x,y) is the transition density of the symmetric a-stable process 
killed on exiting D. 

Lemma 3.1. Let D C R d be an open set which is symmetric with respect to H . 
Then we have 

PD+(t,x,y) = p D (t,x,y) -p D (t,x,y), x,y G D+, t > 0. 
Proof. The proof is based on the reflection principle for the Brownian motion. Put 

t D = mi{s > : (B) Vs £ D}. 

Note that 

r D = inf{s > : B Vs <£ D}, 
f D+ = mi{s > : (B H ) ris <£ D + }, 
Fix x G D + , t > and a Borel set A C D + . We have 

P x (x t eA,f D+ >t) 

= E V P X BH ((B H ) m eA,f D+ >t) 
= E v P£(B Vi eA, Vt <T,T D >t) 

= E V P X B (B vt E A,Tr> > t) — E V P X (B vt G A, rj t > T, r D > t) . (13) 

Note that 

{B vt EA,7] t > T, r D >t} = {B Vt G A, f D > t } . 

Hence ffTBT) equals 

E V P X B (B vt eA,r D >t)- E V P X (B m G A, t d > tj 
= P x {X t eA,r D >t)- P x (Jt t EA,f D > tj 




p D (t,x,y) -p D (t,x,y) dy. 



□ 

Let D C R d be an open set which is symmetric with respect to Hq. If d = 1 < a 
we assume additionally that D is bounded. We define the Green function for D + 
for the process X t by 

POO 

GD+(x,y)= p D+ (t,x,y)dt, x,y G D + , 

Jo 

Go + (x,y) = if x G (D + ) c or y G (D + ) c . For an open bounded set D C R d which 
is symmetric with respect to Hq we define the corresponding Green operator for D + 
by 

G D+ f(x)= G D+ (x,y)f{y)dy. 

J D + 

We assume here that / is a bounded Borel function / : D + — > R. Clearly we have 

G D+ f(x) = E x f D+ f{X s )ds. 
Jo 
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By Lemma 13.11 we obtain the following corollary. 

Corollary 3.2. Let D C R d be an open set which is symmetric with respect to Hq. 
If d = 1 < a we assume additionally that D is bounded. Then we have 

G D+ (x,y) = G D (x, y) - G D (x,y), x,y e D + . (14) 

Lemma 3.3. Let B = B(0, r), r > 0. Assume that f : B — > R is Borel and bounded. 
Then we have 

G B f(x)-G B f(x)= [ G B+ (x,y)(f(y)-f(y))dy. 

J B+ 

Proof. Note that G B {x,y) = G B (x,y) and G B (x,y) = G B (x,y) for any x,y E B + . 
We have 

G B f(x) = / G B (x,y)f(y)dy+ / G B (x,y)f(y)dy 

J B+ J 

G B (x,y)f(y)dy+ / G B {x,y)f(y)dy 

B+ J B+ 



GB(x,y)f(y)dy + / G B (x,y)f(y)dy. 

B+ J B+ 

Similarly, we have 

G B f(x) = / G B (x,y)f(y)dy+ / G B (x,y)f(y)dy 



G B (x,y)f(y)dy + / G B (x,y)f(y)dy. 

B+ J B + 

Using the above equalities and (|14l) we obtain the assertion of the lemma. □ 

Lemma 3.4. Let d > a and V C W C R. d be open sets symmetric with respect to 
Hq. Then we have 

G v+ (x,y) < G w+ (x,y), x,yeV + . 
Proof. Let A C V + be a Borel bounded set. For any x G V + we have 

/ G v+ (x,y)dy = E x H + 1 A (X S ) ds < E* 1 A (X S ) ds = [ G w+ (x,y)dy. 

J A JO JO J A 

Now the lemma follows from ffT4"j) and continuity of Gd{x, ■) on D \ {x} (for D = 
V,W). □ 

Lemma 3.5. Let d > a e (0, 1]. Fix r > and put B = 5(0, r), B + = B + (0,r). 
Then we have 





x — x\ 




\x - y\ 


d—a 


x-y\ 



< G B (x,y) - G B (x,y) < c— — Arl _ n ^ — — , x,y e B+, 
where c = c(d, a). 



g-HARMONIC FUNCTIONS OF FRACTIONAL SCHRODINGER OPERATOR 



Proof. Using Corollary I3.2[ Lemma 13.41 and (j^J) we obtain 
G B (x,y) - G B (x,y) = G B+ (x,y) 

< G Rd+ (x,y) 
= G Rd (x,y) - G nd (x,y) 

= A(d, a) 



\x 


-y\ 


d 


—a _ 




x — 


y\ 


d—a 




x — 


y\ 


d—a 


X 


■ - y\ 


d 


—a 



One can show that for any p > q > and j3 > 

pf> - q P < (2 V2f3)pP- 1 (p-q) 

(we omit an elementary justification of this inequality). Using this one obtains for 
any x,y G B + 

\x-y\ d ~ a ~ \x-y\ d - a < (2 V {2d - 2a))\x - x\\x - y\ d ' a '\ 

which implies the assertion of the lemma. □ 

Now we prove similar lower bound estimates of G B {x, y) — G B (x,y). These lower 
bound estimates will be needed in the proof of Proposition 11.21 We prove these 
estimates only for x G 5 + (0,r/4) and y belonging to some truncated cone lying 
inside B + (0,r). This will be enough for our purposes. The lower bound estimates 
are based on the results of R. Song [27] . 

Lemma 3.6. Let a G (0, 1], d > a. Fix r > and for any x G -B+(0, r/4) put 

K(r,x) = {y= (y!,...,y d ) e R d : r/2 > y { > 2\x\, \y\ < V2 Vi }. 
For any x G -B + (0,r/4) and y G K(r,x) we have 

G mr) (x,y)-G Bm (x,y) > c - ^_ a ^ _ - , 

where c = c(d, a). 

Proof. Let be the ci-dimensional Brownian motion killed on exiting a connected 
bounded open set D C R d , rjt - a/2-stable subordinator starting from zero inde- 
pendent of B® and put = (B )^. Z® is a Markov process with the generator 
— (— A\[)) a / 2 where A\ B is the Dirichlet Laplacian in D. The process Zf has been 
intensively studied see e.g. [27], [28]. Let G B (x,y) be the Green function of the set 
D of the process '. Note that X t = (B H ) Vt = Zf. It follows that 

G B+ ( ,r)(x,y) > G^ +{0r) (x,y), x,y G B + (0,r). 

Let us first consider the case r — 1. Let us fix an auxiliary set f/cR d such that 
U is an open, bounded, connected set with C 1 ' 1 boundary satisfying B + (0, 9/10) C 
U C B + (0, 1). We need to introduce the auxiliary set U because B + (0, 1) is not a 
C 1 ' 1 domain for d > 2 and results from [27J which we use are for C 1,1 domains. 

By [27J Theorem 4.1] we have for any x,y G t/ 

GB(o,i)(x,y) - G B ^x){x : y) = G B+{0tl) (x,y) (15) 

> G| +(0)1) (x,y) (16) 

> G*(x,y) (17) 



> ( 'M^)My) A1 

|x — y\ 2 I \x 



\d-a ' 
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where c = c(d,a). Assume that x G B + (0, 1/4) and y G K(l,x). We have 

\x — x\ 



1 > 



2\x-y\ 



Sjj(x) = 5 B+ (o,i)(x) =Xi = -\x - x\ 



\x-y\>\x- y\, 

4 >^(W + M) = 



M»)>T>i(W + M)> 



It follows that 



5u(x)5u(y) > |x - g| 
|x — y| 2 ~~ 32|x — y| 
Using this and f[T3l- TT5]1 we obtain for x G -B+(0, 1/4), y G i^(l,x) 

G B( o,i) (x, y) - G B (o,i) 1/) > C | a ._^d- Q ^_ y | ' ( 19 ) 
where c = c(d, a). 

Now let r > be arbitrary. Assume that x G B + (0,r/4) and y G K(r,x). Note 
that x/r G -B+(0, 1/4) and y/r G K(l,x/r). By scaling and ffT5j) we get 

Gb(o,i) (->-)- G B(0,1) 



x y I •* 01 I x y I 









x - y| d ~ a |. 


£-y| 



□ 

To obtain estimates of G B+ (x,y) for d = a = 1 we do not use probabilistic 
methods but we use the explicit formula for the Green function for an interval. 

Lemma 3.7. Let d = a = 1. Fix r > 0. For x G R let x = —x. Put B = (— r, r) 

and B + = (0, r). T/ien /or any x, y G 5+ we /love 

< y) -G B (£,y) < ^min (^L, log (^f)) ■ (20) 

For any x, y G (0, r/2) we have 

G B {x,y) - G B {x,y) > -min f — ^ — r,log^ '' 7 



7r \15|x — y|' \4|x — y| 
For any x G (0, r/4) and y G (2x,r/2) we have 

2 \x\ 

G B {x, y) - G B {x, y) > —- -. 

157T \x — y\ 

Proof. By scaling we have 
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so we may assume that r — 1. We have [6] 

G B {x, y) = - log (y/w(x,y) + ^/l + w(x,y)^j , 

where 

, , (i-|*l 2 )(i-M 2 ) 

w{x > y) = — — • 

Let x,y G B + = (0,1). Put t 2 = w(x,y), t\ = w(x,y). Note that t 2 > t\. It 
follows that 

< G B (x, y) - G B (x, y) = - log ^ v ' 2 v 1 ' 2 



= 1 log (l + ^ + ^ ^ .(21) 

It is elementary to show that a/1 + £2 — VI + ^1 < Vh — Vh- Hence ( |2T|) is bounded 
from the above by 



VT V / *l + VI +*1 ^V^ 

2\x + y\ ( 1 



7r \ \x — y\ \x + y\ 



(23) 



< (24) 

Now assume that x G (0,1/4), y G (2a?, 1/2). We will show that G B (x,y) — 
G B (x, y) > 2 |x| / (157r |x — y\). Note that 4|x|/|x — y\ < 4. It is elementary to show 
that for < z < 4 we have log(l + z) > z/5. Using this and (1221 - 124")) we obtain 
that (I2TI) is bounded from below by 

5tt % /i7 +v /r^ 

Note that (1 - |x| 2 )(l - \y\ 2 ) > 1/2 > \x - y\ 2 /2, so 2t x > 1, which implies 2y/t[ > 
VI + £i- Hence f )25|) is bounded from below by 

1 Vh~ Vh \x + y\ ( 1 1 \ 2|x| 



57r 3V^i 15vr Vl 2 - - v\ \ x + y\J 15n\x — y\ 

Now again let x, y G -B+ = (0, 1). We have 



G B (x, y) - G B {x, y) = - log (" v^+VT+^ ( f2G) 

7T VV^l + Vl+^l, 

< - log 



V VTT^ 



\ X -y\2 + {l- 


X 


2 )(i-M 2 ) 


\x + y\ 2 + (1 - 


X 


2 )(i-M 2 ) 



7T I V |x — y\' 

One can easily show that |x + y| 2 + (1 — |x| 2 )(l — \y\ 2 ) > 1 and \x — y\ 2 + (1 — |^| 2 )(1 • 
\y\ 2 ) < 1. Hence ( 127|) is bounded from above by 

1 {2\x + y\ 
log 



7T \ \x — y\ 
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Now let x,ye (0, 1/2). By flU we obtain 



G B (x, y) - G B (x , y) > - log ^ v J 



7T ° \2T/T+h 




\x-y\ 


2 + (l- 


x\ 


2 )(1- 


y\ 


2 ) 


\x + y\ 


2 + (l- 


x\ 


2 )(1- 


y\ 


2 ) 



(2f 

One can easily show that |x + y| 2 + (1 — |x| 2 )(l — \y\ 2 ) < 2 and \x — y\ 2 + (1 — |^| 2 )(1 ■ 
M 2 ) > 1/2- Hence ( 12 8 p is bounded from below by 

1 ( \x + y\ 
-log 



7r \4|x — y\ 



□ 



The estimates of the Green function obtained in this section are crucial in proving 
the main result of this paper. To get these estimates in the transient case we 
used probabilistic methods. There is alternative way of obtaining these estimates. 
Namely, one can use explicit formulas for the Green function of a ball for symmetric 
ct-stable processes (in fact this formula was used in the case d = a = 1). We 
decided to use probabilistic methods instead of explicit formulas for two reasons. 
First, the probabilistic methods are much simpler. Secondly, it seems that it can 
be generalized to some other processes, which are subordinated Brownian motions. 
Especially interesting in this context is the relativistic process, which generator is 

— (y— A + m 2 — m), see e.g. [25], [23], [13]. This operator is called relativistic 
Hamiltonian and is used in some models of relativistic quantum mechanics see e.g. 
[24] . For the relativistic process the explicit formula for the Green function of a 
ball is not known, but it seems that the probabilistic methods from this paper 
could be used to study Schrodinger equations based on the relativistic Hamiltonian 

— (v— A + m 2 — m). 

4. Proof of the main result 
We will need the following technical lemma. 

Lemma 4.1. Fix r G (0, 1], i G {l,...,d} and z = (z\, . . . , z d ) G R d . For any 

x = (xx, . . . ,Xd) G R d denote x = x — 2ej(xj — Zi). Put B = B(z,r). Assume that 
the function f : B — >■ R is Borel and bounded on B and satisfies 

\f{x) - f{x)\ < A\x-zf, xeB(z,r/2), (29) 

for some constants A > 1 and (3 > 0. 

// a G (0, 1) and G [0, 1 — a) then there exists c = c(d, a, (3) such that for any 
x G B we have 

\G B f(x) - G B f(x)\ < cA\x - zf +a + c SnVyeBlfiy)l \x - zf +a . (30) 

r 

// a G (0, 1] and (3 > 1 — a then there exists c = c(d, a, (3) such that for any x G B 
we have 

\G B f(x) - G B f(x)\ < cA\x -z\+ c SUp ^ B \ f ( y >\ x - z y (31) 
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If a = 1 and /3 — then there exists c = c(d) such that for any x E B we have 

\G B f(x) - G B f\x)\ <cA\x- z\^ + c SUIVbI/Q/)I |x - z\W (32) 

r 

Proof. Put B + = {y = (y x , . . . ,y d ) G B : y 4 > 0}. We may assume that z = and 
x = (xi, . . . , Xd) G By Lemma [3.31 we have 



|G B /(x)-G fl /(x)|< / G s+ (x,y)|/(y)-/(y)|dy, 

where Gs + (x, y) = Gb{x, y) — Gb(x, y). We will consider two cases: case 1: d > a G 
(0, 1], case 2: cf = a = 1. We will often use the fact that r G (0, 1] and |x| < r < 1. 

Case 1: d > a G (0, 1]. 

Note that 

|x — x| < 2\x\. (33) 

For any y G 5 + we have 

|x — y| > |x — y|, |x — y| > \x — x|/2. (34) 

Put U x = B(x, \x\) H {y G B + : |y| < r/2}, U 2 = 5 c (x, \x\) (~) {y £ B + : \y\ < r/2}, 
£7g = {y g 5+ : |y| > r/2}. By Lemma [3751, (J2HJ), d33D, flM} we obtain 

G B+ (x,y)\f(y)-f(y)\dy<c\x-x\ [ , l/(y ?~{fr }l , dy 



\x — y| d_a |x — y| 



0, 





x — X 


dy 


x - 


■y\ d ~ a 


\x - y\ 



= I + 11 + III, 

where c = c(d, a). 

For y G U\ we have |y| < |y — x| + |x| < 2|x|. Hence 



I < cA\xf I , = cA\x\ a+ ^ (35) 



\x-y\ 



where c = c(d, a, f3). When |x| > r/4 we get by 

III < csup|/(y)| f = c sup \f(y)\r a <csnp\f(y)\\x\r a -\ 

y&B JB(x,2r) \ X ~ V\ y&B y£B 

where c = c(d, a). If |x| < r/4 then U 3 C B c (x, r/4) n 5(x, 2r) and by (13311 . (J3JJ we 
get 

III < csup |/(y)||x| / | ^— — < csup \f(y)\\x\r a - 1 , 

yeB J B c (x,r/4)nB(x,2r) \ x ~ V\ y&B 

where c = c(d, a). Recall that r < 1. It follows that for x G B + we have 

III<-sup|/(y)||x|, (36) 



r y eB 

where c = c(d, a). 

For y G U 2 we have |y| < |y — x| + |x| < 2|y — x|. Hence 



dy 
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where c = c(d,a,/3). If a E (0,1), (3 E [0, 1 - a) then II < cA\x\ a+l3 , where 
c = c{d,a,P). This, ([H, Q3B} imply ((30]). If a E (0, 1], > 1 - a then II < 
c ^| x | r a+/3-i < cA \ x ^ where c = c ^,a,fi). This, ([35]), (EHD imply (EJJ) . If a = 1, 
/3 = we have 

II = cA|x| / , — . . < cA\x\ [ p~ l dp < cA|d(|log(2r)| + lloglxll) < cA\x\ 1/2 

Ju 2 \x-y\ J\x\ 

where c = c(d). This, ([35]), ([36]) imply (J32]) . 
Case 2: d = a = 1. 
Subcase 2a: x E (0,r/4). We have 

G B+ (x,y)\f(y)-f(y)\dy 

B+ 

2x i-r/2 



< 



rzx pr ■/ 'Z rr 

A{2xf / G B+ (x,y)dy + A / y p G B+ (x, y) dy + 2(sup |/(y)|) / G B+ {x,y)dy 

JO J2x y&B Jr/2 



= I + 11 + III. 

By we have 

f 2x ( Qx \ 
I < cAx 13 / log I r ) dy < cAx^ix + x\ logxl), 

Jo \\ x -y\J 

where c = c(/3). By ([20]) we also have 

r /2 v^dv r /2 

II < cAx / j- — — < cAx / y p ~ x dy, 

J2x F — 1/1 </2:r 

where c = c(/3). Since a; E (0,r/4) by ([20]) we also get 

III < c(sup \f(y)\)x / 1 1 < c(sup \f(y)\)x, 

y&B Jr/2 \ x - U\ V&B 

where c is an absolute constant. 
Now, if > then 

I + 11 + III < cAx + c- sup \f{y)\, 

r yeB 

for some c = c(/3). If P — then 

1/2 

I + II + III < cAr 1/2 + c sup I, 

r v eB 

where c is an absolute constant. 

Subcase 2b: x E (r/4,r). We have 

G B+ {x,y)\f{y)-f{y)\dy 

B+ 

rr/2 pr 

< cAx 13 \ G B+ (x,y)dy + c(sup|/(y)|) / G B+ (x,y)dy 

JO y&B Jr/2 

= I + 11, 

where c = c(P). 

By ([20]) for any x,y E B + we have 

G B+ (x,y) < c-loglar-yl, (37) 
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where c is an absolute constant. 
By fl37|) we obtain 

I < cAx p 



cr 



r/2 



log \x — y\ dy 



( f r/2 \ 

< cAx 13 cr — 2 / log y dy ) 

= cAx^(cr + r — r log(r/2)) 

< cAr^ (r + r | log r | ) , 
where c = c(/3). Similarly, by ( 137|) we obtain 

II < c(sup|/(j/)|) (cr- / log |rr- j/l ) <c(sup|/(y)|)(r + r|logr|), 

y&B V ./r/2 / yes 

where c = c(/3). 
If > then 



I + II < cAr + csup < cAx + c— sup \f(y)\, 

y£B r y£B 



for some c = c(/3). If (3 = then 



I + II < cA(r + r| logr|) + csup \f(y)\ < cAx 1 ^ 2 + c- — sup \f{y)\, 
where c is an absolute constant. 



1/2 



□ 



Lemma 4.2. Fzx r G (0, 1], i G {1, . . . , d} and z = (zi, . . . , zi) G R d . For any 

x = (xi, . . . ,Xd) G R d denote x = x — le-iixi — zA. Put B = B(z,r). Assume that a 
Borel function f satisfies 

\f{x)-f{y)\<A\x-y\\ x,y G B, 

for some A > and T) G (1 — a, 1] . If d > a G (0, 1] £/ien /or any e G (0, r] we have 



B(z,e) 



-^-G B (z,y) 

OZi 



\f(y)-f(y)\dy<cAe^ + 



a-l 



for some c = c(d, a,r]). If d = a = 1 then for any e G (0, r] we have 



\f(y) - f(y)\dy <cAe«(l + \loge\), 



'B(z,e) 

for some c = c{rj). 

Proof. By [TQ1 Corollary 3.3] we have 



d_ 

dzi 



Gb(z,v) 



\ z ~y\ 



\z — y\ A r 

By the assumption on / we have for y G B 

\f(y)-m\<Ay-y\ v = 2 n A\yi-Zi\ v 

If d > a G (0, 1] for any y G B we have 
9 



y G fi, y^z. 



dz. 



G B (z,y) 



\m - m\ < ca ^ b{z jI <ca\ z - yr*- 1 -*, 



\ z - y\ 



1/(2/) - f(v)\ < cA G ^ B _ {z j} v <cA\z- yr 1 log(l + \z - y^ 1 ) 
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for some c = c(d, a, r/). 

If d — a — 1 we obtain from p, Corollary 3.2] that for any y E B we have 

d 

for some c = c(r]). The above estimates imply the assertion of the lemma. □ 

Lemma 4.3. Let a G (0, 1]. Fixr G (0, 1], z = (zi, . . . ,z d ) G R d and i G {1, . . . , d}. 

Put B = B(z, r). Assume that f is bounded and Holder continuous in B with Holder 
exponent r\ G (1 — a, 1], that is 

\f{x)-f{y)\<A\x-y\\ x,y G B. 

Then VG B f(z) exists and we have 

0^G B f(z) = J b q^Gb(z, y)(f(y) - f(y)) dy, (38) 

where B + = {(y u . . . , y d ) G B : yi - Zi > 0} ; and y = y - 2(y, - z^d for y = 
(yi, . . . , yd) ■ We also have 

\VG B f{z)\ < cAr^ +a -\l + | logr|), (39) 

where c = c(d, a, r/). 

Proof. Let g(y) = f(y) — f(z). By our assumption on / we obtain 

\g{y)\<A\y-z\\ y G B(z, r). (40) 
Let h G (-r/8,r/8). We have 

G B f(z + ei h)-G B f(z) = (G B l B {z + eih)-G B l B (z))f(z) 

+G B g(z + eih) - G B g(z). 
By a well known [Tj5] explicit formula for G B 1 B (x) we get 

lim ~(G B l B {z + e t h) - G B l B {z))f{z) = f{z)^-G B l B (z) = 0. 
/i-H) n ozi 

We also have 

I If 

-(G B g(z + dh) - G B g(z)) = - {G B (z + efay) - G B (z,y))g(y) dy 

II n J B{z,2\h\) 

1 

+ h 



(G B (z + e { h, y) - G B (z, y))g(y) dy 

B(z,r)\B(z,2\h\) 



= I + 11. 

We will consider 2 cases: 1: d > a, 2: d = a = 1. 
Case 1: d > a. 

By (140]) and the standard estimate G B (x,y) < K a (x — y) we obtain 



< A2 r >\h\ r >- 1 [ G B {z + e i h,y) + G B {z,y)dy 

JB(z,2\h\) 

< cA^' 1 [ \z + e { h - y\ a - d + \z - y\ a ~ d dy 

J B(z,2\h\) 



< cAlh]^ - 1 , 

where c = c(d, a, if). By our assumption on 77 it follows that lim^o 1 = 0. 



g-HARMONIC FUNCTIONS OF FRACTIONAL SCHRODINGER OPERATOR 
We also have 



17 



II 



dG 



B 



dzi 



z + eih6,y)g(y) dy, 



(41) 



' B(z,r)\B(z,2\h\) 

where 9 = 9(y, z, h, i, a, d, r) G (0, 1). Note that for y G B(z, r) \ B(z, 2\h\) we have 
\y — (z + eih6)\ > \y — z\/2. Using this, (jlDjl and [TQl Corollary 3.3] we obtain for 
y G B(z, r) \ B(z, 2\h\) and 9 as in flU) 



dG, 



dz. 



(z + eih9,y)g{y) 



< cA\y - z\ a+r '- d -\ 



where c = c(a,d,rj). Note that by our assumption on rj the function y — > \y — 
z \a+r]-d-i j g j n tegrable on B = B(z,r). By the bounded convergence theorem we 



get 



It follows that 



lim II 



d 

—G B (z,y)g(y)dy. 

B OZi 







d 

—G B f(z)= I —G B (z,y)g(y)dy 



dzi 



B 



dz. 



(42) 



Note that £-G B {z,y) = -£-G B (z,y), yeB. This and @2) implies ([38 



Case 2: d = a = 1. 
Recall that h G (— r/8, r/8) and r G (0, 1]. By [9j Corollary 3.2] we have 

G B (x,y) < c(l + |log|x-j/||), x,y e B,x ^ y, 



(43) 



where c is an absolute constant. 
By (jUJ) we obtain 



III < ^"l/il"- 1 



G B (z + /i, y) + C7b(2, y) 



B{z,2\h\) 



< cA|/i|' 7 1 / 1 + | log 1 2 + h — y\ | + | log \z — y\ \ dy 

JB(z,2\h\) 

< cA|/i|"(l + |log|/i||), 

where c = c(rf). By our assumption on rj it follows that lim/^o 1 = 0. 
We also have 

C dG 
11=/ —^-(z + h9,y)g(y)dy, (44) 

J B(z,r)\B(z,2\h\) az 

where 9 = 9(y, z,h,r) G (0,1). Note that for y G B(z,r) \ B(z,2\h\) we have 
\y-(z + h9)\ >\y- z\/2. Using this, (HO}. (H5j> and [lT)| Corollary 3.3] we obtain 
for y G r) \ B(z, 2\h\) and ^ as in (jH 



dG } 



dz 



[z + h9,y)g(y) 



< cA\y - z\ r > + \ log \y - z\ 



where c is an absolute constant. Note that by our assumption on rj the function 
y — > \y — z^^il + | log \y — z\\) is integrable on B = B(z,r). By the bounded 
convergence theorem we get 

d -G B (z,y)g(y) dy. 



lim II 

fe-fQ 



dz 
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It follows that 



j-G B f{z) = j B ^G B (z, y)g{y) dy. (45) 



Note that f z G B {z,y) = -£G B (z,y), y G B. This and (JMD implies (EHJ). 

This finishes the justification of f l38]) in both cases. Inequality fl39|) follows from 
(ESD and Lemma [O □ 



Lemma 4.4. Lei a G (0, 2) and .D 6e an open set in R d . For every function f which 
is a-harmonic in D we have 



\Vf(x)\<dff^ xeD. 

b D (x) 



The proof of Lemma fl~4l is almost the same as the proof of Lemma 3.2 in [10] and 
is omitted. 

proof of Theorem \l.l\ Fix arbitrary z = (zi, . . . , Zd) G D and i G {1, . . . , d}. Simi- 
larly like in Lemma I4TT1 for any x = (x±, . . . , Xd) G R d put x = x — 2ei(xi — zi). Using 
[9j Lemma 3.5] let us choose r = r (d, a, ||g||oo) G (0, 1] such that for any r G (0, r ] 
and any ball od radius r contained in D the conditional gauge function for that ball 
is bounded from below by 1/2 and from above by 2. 
Let r = (S D (z) Ar )/2 and B = B(z,r). By © we get 

u{x) = f(x) + G B (qu){x), x G B, (46) 

where f(x) = E x u(X TB ). The function / is a-harmonic on B. 

When u is nonnegative on R d by our choice of ro and by (2.15) in [9] we obtain 
f(x) < 2u(x), x G B. By pH Lemma 3.2] it follows that 

|V/(x)| < dpfi < 4d$L < ej^- } x G B(z,r/2), 
o B [x) d B (z) b D {z) A 1 

where c = c(d, a, ||g||oo)- 

If u is not nonnegative on R d but UnHoo < oo by Lemma [4.41 we get 



|V/(s)| < dffZ < c 7 °° x G 2J(*,r/2), 

where c = c(d, a, ||g||oo)- 

Let K = B when n is nonnegative and K = R d when n is not nonnegative in R d 
and Unlloo < oo. It follows that for any x G B(z,r/2) we have 

\E x u(X TS ) - E*u(X TB )\ < c^f^lx - x\ < c=ff^i|. - 4 (47) 

0d(2) A 1 <)d(zJ A 1 

where c = c(d, a, ||?||oo)- 

Let us consider the following inequality 

\u(x) - u(x)\ < c su to J^(l/)l i _ xeB (z,r/2), (48) 

ivy v ;i - A 1 1 v 7 y v 7 

for some /3 G [0, 1] and c = c(d, a, /3, g, n). 
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Note that \x — x\ < 2\x — z\. Recall that r < 1/2 so \x — z\ < 1/2 for x G B(z, r). 
If (1481) holds then for any x G B(z,r/2) we have 

\q(x)u(x) — q(x)u(x)\ < \u(x)\\q(x) — q(x)\ + \q(x)\\u(x) — u(x)\ 

< csup — x\ v + csup \q(y)\ — T^T~r~T~\ x ~ 

y€K yeD Od[Z) A 1 

S ° &( Z )M |X Zl ' 1 yj 

where c = c(d, a, (3, q, r/). Note that if [3 < 1 — a then (3 A r] = (3 and if /3 > 1 — a 
then /3 A r\ > 1 — a. 

Assume now that (|4"5|) holds for some (a G (0, 1), (3 € [0, 1 — a)) or (a G (0, 1], 
/3 G (1 — a, 1)) and c = c(<i, a, /?, g, r/). 

If a G (0, 1), (3 G [0, I - a) then by 012]) and Lemma EE] we obtain for x G B 

\G B (qu)(x) - G B {qu){x)\ < c ^f'^f' lx - z\^, 

d D (z) A 1 



where c = c(d, a, /3, q, r/). If a G (0, 1], /3 G (1 — a, 1) then by fT4^]) and Lemma I4TT1 
we obtain for x £ B 

|G, WW -G, WW |< JBS^ |,-^, 

d D (-2j A 1 

where c = c(<i, a, /?, g, 77). 

Joining this with ( |4"T]) we obtain in view of (PS)]) that if fPIH]) holds for some (a G 
(0, 1), /3 G [0, 1 — a)) or (a G (0, 1], (3 G (1 - a, 1)) and c = c(d, a, (3, q, 77) then 

Kz) - «(£)] < c su to J M ^)I \ x _ Z \(0+«)M x eB(z,r/2), (50) 

for c = c(d, a, /3, g, rf). 

Note that (T48]) holds trivially for [3 = 0. Assume first that a G (0, 1) and ka ^ 
1 — a for any k G IN. Then repeating the above procedure we obtain that (JIB]) holds 
for [3 = 0, a, 2a, . . . and finally for (3 = 1. 

Assume now that a G (0, 1) and k^a = 1 — a for some ko G IN. Then we obtain 
that flU]) holds for f3 = 0, a, 2a, . . . , k a. Then (02]) holds for any f3 G [0,k a\. 
In particular, it holds for (3 = k a — a/2. By floD]) we obtain that (T4"5]) holds for 
(3 = k a + a/2 G (1 — a, 1). Then, again by flBU]) we obtain that PS]) holds for (3 = 1. 

Finally assume that a = 1. (E2J) gives that (ggj) holds for /3 = 1/2. Then by ((SD]) 
we obtain that (J35J holds for /3 = 1. 

Now let us fix arbitrary w E D and put s = (5d(w) A r )/8. We will show 
that Vu(w) exists. Let us take x,y G B(w,s). Since z £ D was arbitrary one 
can take z = (x + y)/2 and choose the Cartesian coordinate system and i so that 
y = x = x — 2ti{xi — zi). We put r = (<5d(<2) A r )/2 as before. Note that 

or>U) > od(w) - — = — - — > 7s, s < —. 



We also have 

I I \ x ~y\ ^ / ( <W) r \ r 
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so x G B(z, r/2). On the other hand we have Sr>(z) < s + Sr>(w) so 

r = 5d{z) A r ° < {S + A r ° < £ + A r ° = -, and \w - z\ < s. 

2 2 ~2 2 2' ii- 

Hence B(z,r) C B(tu, lls/2) which gives sup pgB(zr) |«(p)| < sup peB{w Us/2) \u(p)\. 
By (gHD for p = 1 we obtain 

i / \ f x, . sup |«(p)| sup Kp)|, . 

- M !/ < c X - z < c x - y\, 

do{z) A 1 Od(w) A 1 

where c = c(d,a,q,rj) and iT' = -B(if, lls/2) when it is nonnegative in R d and 
K' = R d when w is not nonnegative in R d and ||m||oo < oo. Since x, y G B(w,s) 
were arbitrary we obtain that gw is Holder continuous with Holder exponent rj A 1 
in -B(u>, s). 

Using (P) for W = B(w,s) and Lemma [O] for B(w,s) we obtain that Vm(w) 
exists. Since w G -D was arbitrary this implies that Vw is well defined on D. 

Now again let us fix arbitrary z G I?, i G {1, . . . , d} and put r = (#0(2) A r )/2, 
B = B(z,r), K = B when m is nonnegative and X = R d when u is not nonnegative 
and |H|oo < 00. 

When u is nonnegative, by the Harnack principle (see [3 Theorem 4.1]) we have 

sup \u(y)\ = swpu(y) < cu(z). (51) 

By the proof of [9j Theorem 4.1] it follows that c = c(d, a, ||g||oo)- 
Put x = z + he u h G (0, r/2). By @BJ for /3 = 1 we get 

\u(z + hei) - u(z - hei)\ < c ^ v& ,\ ■ ^- ■ fe, 

^(2) A 1 

where c = c(d, a,q,rj). Since z G {1, . . . , d} is arbitrary it follows that 

ivy / M ^ su P ^|m(|/)| 

|V ^ }I ^ C m*)ai • (52) 

Of course this gives (j3J. When it is nonnegative ( 1521) and (15T|) imply Q. □ 

5. Proof of Proposition 11.21 and Theorem 11.31 

First we prove Proposition 11.21 By saying that Vu(x) exists we understand that 
for each i G {1, . . . , d} \im.h->o(u(x + hei) —u(x))/h exists and is finite. We say that 
a function is Holder continuous if it is bounded and measureable. 

proof of Proposition IT7B . Let us choose arbitrary point w G D and r G (0, <5d(u>)/3). 
Put 

q( x ) = l B ( w , r )(x){r 2 -\x- w\ 2 y- a , x G R d . (53) 

It may be easily shown that q(x) is (1 — a) Holder continuous. We may assume that r 
is sufficiently small so that (D, q) is gaugeable. (The fact that (D, q) is gaugeable for 
small r follows by Khasminski's lemma, see page 57 in [8].) Put u(x) = E x (e g (ro)) , 
x G R d . u(x) is the gauge function for (D,q). By Theorem 4.1 in [9] u(x) is regular 
g-harmonic in D. Note that u is continuous and bounded on D. 

Fix z G dB(w, r). We may assume that the Cartesian coordinate system (xi, . . . , 2^) 
is chosen so that z = (0, . . . , 0) and w = (r, 0, . . . , 0). Let B = B(0, r). We will 
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show that Vw(0) does not exist. On the contrary assume that Vw(0) exists. By ([9]) 
we have 

u{x) = E x (u(X TB )) + G B (qu)(x) } xeB. (54) 

Of course V ' E x {u(X TB )) exists for x G B (see (10) and Lemma 3.2 in [ID]). Put 
fo(y) = u(0)q(y) and f x {y) = (u(y) - u(0))q(y). We have 

u(y)q(y) = fo(y) + fi(y). (55) 

For any s > put B+(0,s) = {( yi ,...,y d ) G B(0,s) : y x > 0}, B + = B + (0,r) 
and y = y — 2e x y x for y = (y x , . . . , y^) G R d . Recall that we have assumed that 
Vtt(O) exists. By this and boundedness of u we get \u(y) — u(0)\ < c\y\ for some 
c = c(w, z, r, d, a, q) and any y G R d . It follows that for y G -B+(0, r/2) 

l/i(2/) - = \fi(y)\ = \u(y) - u(0)\\q(y)\ < c\\q\Uy\, 

where c = c(w, z, r, D, d, a, q). By (l3Ti) for any x G 5(0, r/2) we have 

|G fl /!(x) - G B /i(x)| < c\x\ + c SUP?/g ^ l/l(y)l |x| < c\x\, (56) 

for some c = c(w, z, r, D, d, a, q). 

Now let us consider the case a G (0, 1], d > a. By Lemmas 13 . 3 1 and [3761 we get for 
xeB + (0,r/4) 



G B fo(x) - G B fo(x) = / (G B (x,y)-G B (x,y))(f (y)- f Q (y))dy 



= u(0) / (G B (x,y)-G B (x,y))q{y)dy 
Jb + 

> CM (°) / .., g(y)^, (57) 



where c = c(u>, z, r, D, d, a, q) and K{r, x) is defined in Lemma 13.61 for i = 1. Since 
it is positive on D we have u(0) > 0. 

Note that for x G _B + (0,r/4) and y G if(r,x) we have \x — y\ < (3/2)|y|, \x — y\ < 
(3/2) I y I . Hence f )57|) is bounded from below by 



cu{0)\x-x\ / \y\ a d 1 q(y)dy, 

JK(r,x) 

where c = c(w, z, r, D, d, a, q). One can easily show that for any x G B + (0, r/4) and 
y G K(r,x) we have q(y) > c|t/| 1_q , where c = c(d, a, r). Let x = (x\,0 . . . ,0) G 
B+(0,r/4). It follows that 

G B f (x) - G B f (x) 



>cu(0) \y\- d dy, (5* 



where c = c(u>, r, D, d, a, q). It is clear that if x = (x%, . . . , 0) tends to then the 
right-hand side of (155]) tends to 00. This, (|54|) . the fact that VE x (u(X TB )) exists for 
x G B, (155]) and ( 156]) give contradiction with the assumption that Vtt(0) exists. 
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Now we will consider the case d — a — 1. Recall that in this case q(y) 
1-B(w,r)(y) = l(o,2r)(|/)- By Lemmas 13.31 and [3771 we get for x G (0, r/4) 

G B f (x)-G B f (x) u(0) 



\x — x\ 2\x\ 



(G B (x, y) - G B (x, y))q(y) dy 



> u(0) f r/2 dy 
15vr J 2x y ' 

It is clear that if x — > then (15$]) tends to oo. This, (151]) . the fact that -^-E x {u(X TB )) 
exists for x £ B, ( 155]) and ( 156]) give contradiction with the assumption that u'(0) 
exists. □ 

Now we will prove lower bound gradient estimates. The idea of the proof is to 
some extent similar to the proof of lower bound gradient estimates in [10]. The 
main difference is the use of Lemma 15.11 below instead of [10] Lemma 5.4]. There 
are essential differences in proofs of Lemma 15.11 and [TO] Lemma 5.4]. The key 
arguments in the proof of Lemma 15.11 are based on Lemma 14.31 

We will use the notation as in [TU]. For x = (xi, . . . , Xd) G R d we write x = (x, Xd), 
where x = (x\, . . . ,Xd~i)- In order to include the case d = 1 in the considerations 
below we make the convention that for x G R, x = and we set R° = {0}. 

We fix a Lipschitz function T : R d_1 — > R with a Lipschitz constant A, so that 
|r(x) — T(y)\ < X\x — y\ for x,y e R d_1 . We put p(x) = Xd — r(x)- D denotes the 
special Lipschitz domain defined by D = {x 6 R d : p(x) > 0}. The function p(x) 
serves as vertical distance from x G D to dD. We define the "box" 

A(x, a, r) = {y G R d : < p(y) < a, \x — y\ < r}, 

where x G R d and a, r > 0. We note that A(sc, a, r) is a Lipschitz domain. We also 
define the "inverted box" 

V(x, a, r) = {y G R d : -a < < 0, |x — y\ < r}. 

The same symbol V is used for the gradient but the meaning will be clear from the 
context. 

For r > and Q G dD we set A r = A(Q, r, r) and G r = G& r . For a nonnegative 
function u we put 

u Ar (x) = E x u(X TAr ) } x G H d . 
Fix Q G and assume that a Borel function q satisfies 

\q{x)-q(y)\<A\x-y\\ (60) 

for some A > 0, rj G (1 — a, 1] and all x,y £ A So for some so G (0, 1]. 

Now we will repeat the assertion of Lemma 5.3 [10]. Note that the assertion of 
Lemma 5.3 in [TU] holds for all a G (0, 2) under the condition that ^Ia^ G £[ q for 
some So G (0,1]. This condition follows from ( I60p . 

For every e > there exists a constant r = r (d, \,a,r],q, s ,e) < s < 1 such 
that if r G (0, r ] and u : R d [0, oo) is g-harmonic and bounded in A r then 

(1 -e)u Ar (x) < u(x) < (1 + s)u Ar (x), x G R d , (61) 

and 

G r (|g|u)(aO < £M Ar (x), x G R d , (62) 
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Lemma 5.1. Let a G (0, 1] and e G (0, 1/2]. There exist constants c = c(d,a,rj,q) 
and k = n(d, A, a, r], q, r , Sq, e) < r swc/i t/iai if < r < k, u : R d — > [0, oo) is 
q-harmonic and bounded in A r then 

\VG r (qu)(x)\ < ec x G A r . 

SaAx) 



Proof. Let us choose 



k = max{s G (0, r ] : sup a'' (1 + | loga|) < e}. 

0<a<s 

Fix r G (0,/c] and x G A r . Note that 5A r {x ) < r < 1. Let B = B(a;o, 5& r {x Q )/2). 
We have 

G r (gu)(x ) = G B (g M )(a;o) + £r°G r (<H(*r B ). (63) 

We will estimate gradient of two terms on the right-hand side of ( 16 3 p separately. 

Let Pb{x, z), x G B, z G int(i? c ) be the Poisson kernel for B (that is the density 
of the P x distribution of X{r B ) [S])- By Lemma 3.1 in [TU] we have 



|V (E x °G r (qu)(X. 



1~B> 



V z / P B (x , z)G r (qu)(z) dz 



B c 



< / \V x P B (x ,z)\G r {\q\u){z)dz 



c 



< t — t — r / P B (x ,z)G r (\q\u)(z) dz 

OaAXo) 



B c 



= ^-rE*°G r (\q\u)(X TB ), (64) 

where c = c(d, a). 

By fl63|) for \q\ instead of q we obtain that E x ° G r (\q\u)(X TB ) < G r (\q\u)(x ). Using 
this and f lBTj) . flB2J) we obtain that is bounded from above by ceu{xo)/5^ r {xo), 
where c = c(<i, a). 

By Theorem II. II for any x, y G B we have 

|m(x) - < 7 r ( supw(z) J |x - y\, 
Oa t {x ) \zeB J 

for some c = c(d, a, rj, q). Using this and flBUl) for any x,y £ B we get 

|g(x)w(x) - g(y)u(y)| < \q(x)\\u(x) - u{y)\ + \q(x) - q(y)\\u(y)\ 

c(A + l) ( . \ , 
< — — - t — — ( sup u(z) I \x — y\ ', 
zeB J 



SaAxq) 



for some c = c(d, a, rj, q). 

Hence by Lemma [4.31 we obtain 



\VG B {qu){x Q )\ < C [ A i 1 } ( supu(z) ) (^(xo))^-^! + | Iog(«J Ar (xo))|), (65) 



^A r (a;o) 



where c = c(d, a,r],q). 

Note that SA r {xo) < r < k. By our choice of n we have 
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Using this and the Harnack inequality (see (j5T]) with y changed to z and z 
changed to xq) we obtain that the right-hand side of f[65l) is bounded from above by 
ceu(xo)/SA r (xo), where c = c(d,a,r],q). □ 

The next lemma is similar to Lemma 5.6 in [TO] . 

Lemma 5.2. Let a £ (0,1]. There are constants c = c(d,a,X), h = h(d,a,X) 
and ri = ri(d,a, \,rj,q, sq) such that if < r < r\ and u is nonnegative in H d , 
q-harmonic and bounded in A r and vanishes in V(Q,r, r) then 

u( x) 

\Vu(x)\ > c r y . , x £ A(Q,rh,r/2). 
5 Ar {x) 

Proof. The function u satisfies ([9]) with W = A r . Using pjjj Lemma 4.5] and scaling, 
(1611) and Lemma 15.11 we obtain the result by an appropriate choice of e in Lemma 



o □ 

Proof of Theorem \1.3[ The upper bound follows from Theorem 11.11 The lower 

bound follows from Lemma 15.21 and compactness of dD D K . □ 

6. Applications 



As an application of the main results of this paper we obtain gradient estimates 
of eigenfunctions of the fractional Schrodinger operator. 

Corollary 6.1. Assume that a £ (0,2), D C R d is an open bounded set, q £ J a_1 
when a £ (1,2), or q is Holder continuous on D with Holder exponent r\ > 1 — a 
when a £ (0,1]. Let {<£> n }^i be the eigenfunctions of the eigenvalue problem ( TJ^ )- 
< T77p for the fractional Schrodinger operator on D with zero exterior condition. Then 
V(p n (x) exist for any n £ IN, x £ D and we have 

\V Vl {x)\<c / f ^l v xeD, (66) 
5 D (x) A 1 

where c = c(D, q, a, rj) and 

|v ^» £ sMXi- x e B ' (67) 

where c n = c n (D,q,a,r]). Furthermore, if additionally D C R d is a bounded Lips- 
chitz domain then there exists e = e(D,q,a,i]) such that 

|V?i(aO| >cf44> 5d(x)<£, (68) 

where c = c(D, q, a, rf) . 

The result is new even for q = 0. In that case this is the eigenvalue problem for 
the fractional Laplacian with zero exterior condition. This eigenvalue problem have 
been recently very intensively studied see e.g. [JJ, [U], [2], [21], [18], [3]. 

For a = 2, under additional assumptions that <i > 3, D is connected and Lipschitz, 
inequalities (1661) . (167|) follows from [TOJ Theorem 1] and inequality (168]) follows from 
[U Theorem 1]. 

Before we come to the proof of Corollary 16.11 we will need the following easy 
addendum to the results obtained in [TO] . 
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Lemma 6.2. Let q 6 (1,2), g £ d"^ 1 and D C R d be an open set. Assume that the 
function u is q-harmonic in D and ||w||oo < oo. Then V«(i) exists for any x G D 
and we have 

\\u\\ 

\Vu{x)\ < c ' °° - , xED, 

where c = c(d, a, q). 

Proof. The proof of this lemma follows from the arguments used in [TO] . First note 
that the assertion of Lemma 5.4 in [10] remains true if we replace the assumption 
that u is nonnegative in R d by the assumption that ||w||oo < oo and when we replace 
u(x) by ||w||oo on the right-hand side of the estimate of \VG r (qu)(x)\. Then the 
proof of Lemma [6.21 is almost the same as the proof of Lemma 5.5 in [TU]. □ 

proof of Corollaru \6.1\ It is clear that (p n is not (q + A n )-harmonic on the whole D 
because (D, q + \ n ) is not gaugeable. However by the definition of the Kato class and 
standard arguments (see e.g. page 299 [9]) for any n — 1, 2, . . . there exists r G (0, 1] 
and the finite number of balls B(x\, r), . . . , B(xm, r) such that X\, . . . , xm G D, 

M 

D C J2 B ( x m,r) 

m=l 

and each (B(x m , 2r) fl D, q + A n ) is gaugeable. This means that <p n is (q + A n )- 
harmonic on each B(x m , 2r) fl D. Note that for any x G B(x m , r) fl D we have 

8 B (x m ,%r)nD{x) A 1 > 8 D {x) A r A 1 > r(S D (x) A 1). 

Now, fl66|) . fl67|) follow from Theorem 11.11 for a G (0,1] and from [TOj Lemma 5.5], 
Lemma 16721 for a G (1, 2). Inequality fl68|) follows from similar arguments and Lemma 
ETJfor a G (0, 1] and (TUl Lemma 5.6] for a G (1, 2). □ 

As another application of our main result we show that under some assumptions 
on q a weak solution of A a ^ 2 u + qu = is in fact a strong solution. First we need 
the following easy lemma. 

Lemma 6.3. Let a G (0,1). Choose x$ G R d and r > 0. Assume that a Borel 
function u : R d — > R satisfies 

f Hy)\ 

/ 1 , , ay < oo, 

Vm(x) exists and \Vu(x)\ < A for all x G B(xo,r) and some constant A. Then 
A a / 2 u(x) is well defined and continuous on B(xo,r /2) . 

Proof. Let x G B(x ,r/2). Choose e G (0,r/2). We have 

\u(y) — u(x)\ , T Al-y — x\ , . . 

- \ d J d y ^ / ,, (69) 



j/|<£ 12/ X \ J\x-y\<e \V 



cA I p~ a dp 0, when e ->• 0, (70) 
Jo 



where c = c(d). By the definition of A a ^ 2 (see Preliminaries) we obtain that 
A a / 2 u(x) is well defined. 

One can easily show that for any fixed e G (0, r/2) the function 

/.<*>=/ 



X — J/|>£ 



y - s 
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is continuous on B(xo,r/2). This and ( 169|) - (!7D|) imply that A a ^ 2 u(x) is continuous 
on %,r/2). □ 



proof of Corollary 1.4 Choose arbitrary xq G -D. It is clear that there exists r > 
such that B(xo,2r) CC D and (£>(a;o, 2r), g) is gaugeable. This can be done by 
Khasminski's lemma (see page 299 [9]). Put B = B(xo,r). By Theorem 5.5] we 
may assume that u is a g-harmonic function on B(xo, 2r) (after a modification on a 
set of Lebesgue measure zero). By ([9]) we get 

u(x) = E x u(X TB ) + G B (qu)(x), x G R d . (71) 

By Theorem 1 1 . 1 1 and Lemma l63l A a//2 M(x) is well defined and continuous on B(xq, t/2). 
The function v(x) = E x u(X TB ) is an a-harmonic function on B(xq, r/2), so A a ^ 2 v (x) = 
on B(x ,r/2). Hence by fjTTj) we obtain 

A a / 2 u(x) = A a / 2 (G B (qu)){x), x G B(x ,r/2). 

By Lemma 5.3 [9] we have 

A a / 2 (G B (qu))(x) = -q(x)u(x), 

for almost all x G B(x ,r/2). But both sides of this equality are continuous so in 
fact this equality holds for all x G B(x Q ,r/2). □ 
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